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^ . Abstract 

O I We derive generic formulas for the polarization density matrix of leptons produced 

in uN and VN collisions and briefly consider some important particular cases. Next we 
employ the general formalism in order to include the final lepton mass and spin into the 
^ , popular model by Rein and Sehgal for single pion neutrinoproduction. 
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Polarization of leptons generated in uN and VN collisions is important for studying neutrino 
oscillations and relevant phenomena in experiments with atmospheric and accelerator neutrino 
beams. Let us shortly touch upon a few illustrative examples. 
. • Contained r lepton events provide the primary signature for z/^ Ur oscillations. Be- 

O \ sides they are a source of unavoidable background to the future proton decay experiments. But 

\ a low or intermediate energy r lepton generated inside a water Cherenkov detector is unobserv- 

able in itself and may only be identified through the r decay secondaries whose momentum 
configuration is determined by the r lepton helicity. 
^ . • In case of — mixing, the leptonic decays of r's generated inside the Earth yield 

an extra contribution into the flux of through-going upward-going muons (TUM) and stopping 
\ muons (SM) which is absent in case of Vy^ — Ug or — mixing. The absolute value and en- 

^ ■ ergy spectrum of the "r^s muons" are affected by the r beam polarization. The contribution is 

evidently small but measurable in future large-scale experiments, particularly those with mag- 
netized tracking calorimeters (like in the experiments NuMI-MINOS and MONOLITH). Note 
that the energy and angular distributions of the charge ratio for the "r^3 muons" are consider- 
ably differ from those for the "direct" TUM and SM since the longitudinal polarizations of r+ 
and T~ have in average opposite signs. 

• Decay of z/^ or 77^ induced muons with energy below the detection threshold may produce 
detectable electrons whose energy distributions are affected by the muon polarization. Such 
events, being classified as "e-like" (for a water detector) or "showering" (for an iron detector), 
mimic the z/g or z/g induced events. 

In this paper, we derive general formulas for the lepton polarization density matrix by ap- 
plying a covariant method (Sect. 2) and briefly consider their applications to deep inelastic, 
quasielastic and resonance neutrino interactions. We explicitly demonstrate that the perpen- 
dicular and transverse polarizations are dependent of an intrinsically indeterminate phase and 
thus unobservable in contrast with the longitudinal polarization and degree of polarization. In 
Sect. 3 we discuss with some details a generalization of the Rein-Sehgal model for single pion 
neutrinoproduction through baryon resonances which takes into account the final lepton mass 
and spin. 
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2 Polarization density matrix 

The lepton polarization vector V = {Vp, VtjVl) is defined through the polarization density 
matrix p = | (1 + crP) whose matrix elements are given by contracting the leptonic tensor 
L'^y with the spin- averaged hadronic tensor Wap- The leptonic tensor is given by 



Jx (rf)* with = «(fc',s)7" (^^) «(^) for 
J\ (ly) with = I;(/c)7" ^) ^ (^'' 



where k and A;' are the 4-momenta of or Ve and lepton £~ or = e, r), A and A' are the 
lepton helicities, s is the axial 4-vector of the lepton spin. 

It can be shown that the weak leptonic currents and are expressed as 

= Nx [k'^'iks) - s^'ikk') - ie'^'^'^hpk^k's + rruk'^] , (2a) 

= \Nx [k'^'iks) - s"{kk') + ie"^"'^S0k^k's - mek''] . (2b) 

Here is the lepton mass and the normalization constant Nx is expressed in terms of the 
kinematic variables and of two intrinsically indeterminate phases 99+ and <^_: 



^ ^ exp {±i(px) ^ exp {±i(px) / Ej^ XPj 
^ y/(kk')± me(ks) rnt E,{It XcosO)' 

Here By, Ei and d are, respectively, the incident neutrino energy, lepton energy and the scatter- 
ing angle of the lepton in the lab. frame. Pi — yjE^ — m^; the upper (lower) signs are for ue 
{pe). As it follows from Eq. (2a), the components of the neutrino current are 

fx = exp (z^x) V(l-Acos^)E,(E,- AP,), 
j. = ^[AP,k-E.k' + i(kxkO]. 

As it follows from Eqs. (2) and (3), the currents for neutrino and antineutrino are related by 

jr=-A(A)*. (5) 

We use the standard representation for the hadronic tensor (see, e.g., Ref. [1]) 



which includes 6 nucleon structure functions, Wn, whose explicit form is defined by the par- 
ticular subprocess (QE, RES or DIS). Here p and M are the nucleon 4-momentum and mass, 
respectively, q = k — k'istheW boson 4-momentum. By applying Eqs. (1), (2) and (6), we 
obtain 



, _ Lfx'W^p _ m^E^NxNl, ' , 

k=l 



2 
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where rjx = (1 + A)/2 (77+ = 1, ?]_ = 0) and the dimensionless normaUzation factor TZ is given 
by the condition Trp— 1: 

= ilk ('-+^--) 

' E^ — Pi cos ^ 



± 



M 
M 



Wi + 



2M2 

— Pi COS ^ 

2M 



1^4 + 



El + Pe COS 9 



2M 



Wo 



2M2 



rric 



Now we can find the explicit formulas for the elements of the polarization density matrix in 
terms of variables E^, Pi and 9: 

P++ {E,, Pe, 9) = p__ {E,, -Pe, n - 9) ^ ^M^^^ 

{E^, Pe, 9) = p% {E^, Pe, 9) = {X - zy) e^^. 

Here we have introduced the following notation: 



y 
z 



P^ Tir 



(licos^) (wi± 



E.^Pe 



W: 



+ 



1 ^ cos ^ 



2M 

Ee^ Pe I Ee± Pe 



M 



M 



W4 - W5 



and if = (fij^ — (f^. Finally the projections of the lepton polarization vector are given by 



Vp\ _ mism9 / cos 99 smip\ fX\ 
Pt/ T/Vf^ V - sin ^ coscpjiyj' 



(7a) 
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(7b) 



By putting 93 = / the formulas for Vp and Vl exactly coincide with those of Ref. [2] (obtained 
within a noncovariant approach under assumption Wq — 0). 

Several simple conclusions immediately follow from Eqs. (7). First, the perpendicular and 
transverse projections are unobservable quantities in contrast with the longitudinal projection 
of P and the degree of polarization \ V\. Second, supposing that Wq = (as is the case in the 
Standard Model) one can force the polarization vector to lie in the production plane. Third, a 
massless lepton is fully polarized, V — (0, 0, ^1). In particular, at the energies of our interest, 
electron is always fully polarized while in general, this is not the case for muon and r lepton. 

We end this section with a short description of the structure functions relevant to the three 
fundamental subprocesses, DIS, QE and RES. 

2.1 Deep inelastic scattering (DIS) 

In this case the relation between the structure functions W^^^^ {x, Q"^) and measurable quantities 
Fji {x, Q"^) is obtained in a straightforward manner: 

Wr'{x,Q')=F,{x,Q'), Wr''{x,Q')^w-'F^{x,Q'), n = 2,...,6. 

Here = — g^, x = Q'^/2{pq) is the Bjorken scaling variable and w = {pq)/M'^. 
The generally accepted relations between the functions Fi, F2, F4 and F5 are 



^ 2x{l + R)\ x' 



2x\2x ' ' 



where it! is the ratio of longitudinal to transverse cross sections in DIS and x' — j (4M^). 



2.2 Quasielastic scattering (QE) 

The charged hadronic currents describing the QE processes are written [1] 

(P, PVa 1?^, P) = COS OcUp (p) Fa M„(p), 
(n, p'\J~\p, p) = COS 9cUn {p') FaUp{p), 

where 9c is the Cabibbo mixing angle, p' = p + q and the vertex function 

7-1 771 I • ^3 , Qa . ( 771 I Pa ^ Pa , 7^ 1 

Fa = -faFv + tcJap^FM + —Fg + ( 7a^A + ^ Ft + —Fp ] 75, 

is defined through the six, in general complex, form factors Fj (g^). A standard calculation then 
yields 

^) {x, = cos^ Ocw-'un (Q') S{l-x), n = 1, . . . , 6, 
^ We adopt this convention from here on. Therefore, according to Eq. (7a), Vp oc X and Vt oc y. 



4 



where the functions uin are the bihnear combinations of the form factors: 



uJi = + x' {\Fa\' + \Fv + Fuf) , 
u;2^\Fvf + \FAf + x' (|Fm|' + 4 [Frf] , 
cos^ - 2Rt [FX (Fv + Fm)], 

= Re [F* {Fs - |Fm) - ^ (Ft + Fp)] + x' {\ \Fm - Fgf + \Ft + Fp\^) 
-\{1 + x')\Fm? + {1 + \x')\Fs\\ 
a;5 = 2 Re [F| (Fy - x'Fm) - F* {Fa - 2x'Fp)] + U2, 
a;6 = 2 Im [F| {Fy - x'Fm) + F* (Fa - 2x'Fp)] . 

The only difference between this result and that from Ref. [1] is in the relative sign of the terms 
in loq.'^ Assuming all the form factors to be real we have coq — O and thus Vt — 0. 



2.3 Single resonance production (RES) 

Let us now consider the case of single A resonance neutrinoproduction, 

+ n{p) ^ + A+ (A++) , + n{p) ^ + A" (A°) . 

Assuming the isospin symmetry and applying the Wigner-Eckart theorem, the hadronic weak 
current matrix elements are given by [2] 

(A+, p'\J^\n, p) = (A°, p'\J^\p, p) = cosdc^^{p')r,,su{p). 
(A++, p'\J^\p, p) = (A-, p'\Ja\n, p)^V3 coseci^f" {p') r^^u{p). 

Here ?/^" {p') is the Rarita-Schwinger spin- vector for A resonance, u{p) is the Dirac spinor for 
neutron or proton and the vertex tensor F^p is expressed in terms of the eight weak transition 
form factors C^f^^^^ (Q^) [3]: 



F 



a/3 



75 



+ c. 



A9 apq - laqp 

M 



, r^A 9ap {qp') -p'gqp , ^a„ , ^/ 

^4 ^ <^5 9al3 + <^6 



M2 ■ 



Below we will assume that time-reversal invariance holds so that all the form factors are rela- 
tively real. After accounting for the explicit form of a spin 3/2 projection operator and comput- 
ing the proper convolutions we arrive at the following expressions for W^^^: 



W^'' = ^cos' 9cMM^\r,i^{W)\' {VfCjC^ + Ai^CfC^) , n ^ 3, 

j,k=3 to 6 

W^'^ ^2kcos'9cMM^\v^^{W)\' Y1 ^nCjCt and M^r'^ = 0. 

j,k=Z to 6 



(8) 



^According to Llewellyn Smith, the functions = ^5 — 0J2 and are, respectively, the real and imaginary 
parts of a unique function. Our examination does not confirm this property for the general case of nonvanishing 
second-class current induced form factors Fs and Ft. 
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Here k = 2/3 for A+ and A° production or k = 2 for A++ and A production; 



1 

2^ 



FAiW) 



{w^-Miy + ri{w)/4: 



is the Breit-Wigner factor and / a(W^) is the running width of the A resonance which can be 
estimated by assuming the dominance of -S-wave (L — 0) or P-wave (L = 1) A — > Nn 
decay [3,4]: 







2L+1 









where 



M2 + 



2W 



is the pion momentum in the rest frame of the A with invariant mass W — \p'\ — \p-\- q\; Ma 
and are, correspondingly, the central mass and decay width of A. 

The coefficients V^^, and K^j^ are found to be cubic polynomials in invariant dimen- 
sionless variables x, w and parameter C, — M/Ma- Only 70 among the total 144 coefficients 
are nonzero (see appendix). 



3 Single pion production in Rein-Sehgal model 

The Rein-Sehgal (RS) model [5] is one of the most circumstantial and approved phenomenolog- 
ical tools for description of single-pion production through baryon resonances in neutrino and 
antineutrino interactions with nucleons. It is incorporated into essentially all MC neutrino event 
generators developed for both accelerator and astroparticle experiments and is in agreement 
with available experimental data (see, e.g., Ref. [6] for a recent review and further references). 
However the RS model is not directly applicable to the Vr and Vr induced reactions since it 
neglects the final lepton mass. Due to the same reason, the model is not suited for studying 
the lepton polarization. In this section, we describe a simple generalization of the RS model 
based upon the covariant form of the charged leptonic current jx with definite lepton helicity 
discussed in Sect. 2 which allows us to take into account the final lepton mass and spin. 

The charged hadronic current in the RS approach has been derived in terms of the Feynman- 
Kislinger-Ravndal (FKR) relativistic quark model [7] and its explicit form has been written in 
the resonance rest frame (RRF); below we will mark this frame with asterisk (*). In RRF, the 
energy of the incoming neutrino, outgoing lepton, target nucleon and the 3-momentum transfer 
are, respectively, 

E* = {2MEi + g2 _ mj) = [MEe + E, {E^ - cos 9) - mf\ , 

M 

E*^ - (Et - E^) ^ -{M + E, - E,) , 



M , , M /— : 

T7^h\ = T7^^E-^-2E^Pecose + Pl 
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It is convenient to direct the spatial axes of RRF in such a way that p* = (0, 0, — |q*|) and 
k* = k'y* — 0. These conditions lead to the following system of equations: 



|q I, 



(10) 



By using Eqs. (9) and (10) we find the components of the lepton spin 4-vector: 
1 



s* — 

On — 



meW 



mi |q| 
1 



sm.9, 



y 



rui |q| W 



[{E, cos 9 - Pe) {MEe - m\ + E,E() - E^Pi {E, - Pe cos 6)] 



Then, by applying the general equations (4), the components of the leptonic current in RRF 
with the lepton helicity A measured in lab. frame, are expressed as 

Jo* = N^rne^il - XcosO) {M-E,- \P,) , 
= Nxm,^ sine (Pe-XE,), 



jy = iXNxmiE^ sin 9, 

E, 



Nxmt 



q\W 



1 - A cos e) [{E^ + XPe) {M - Ee) + Pe {XE^ + 2E^ cos 9 - Pe)] . 



On the other hand, in the spirit of the RS model, the leptonic current still can be treated as 

the intermediate W boson polarization 4-vector. Therefore it may be decomposed into three 
polarization 4-vectors corresponding to left-handed, right-handed and scalar polarization: 

e2 = -^(0,l,-^,0), 
eS=^(0,-l,-i,0), 



p°' — 



Here the vectors c'l and are the same as in Ref. [5] while e"^^ has been modified to include 
the lepton mass effect: 



(A) 



|q| 



1 



Then the coefficients are explicitly defined through the components in RRF as 
For antineutrino induced reactions one have to take into account relation (5). It yields 

[QUl=[QU)].- h*A)],= K-A)]„. 

Within the extended RS model, the elements of the polarization density matrix may be 
written as the superpositions of the partial cross sections^ a^^' , a^' and cr^^': 



Paa' = 



-, ^AA' - 

i=L,R,S 



,A^A' AA' 



and the differential cross section is given by 



dQ^dW^ 27r2M|q|' 

The partial cross sections are found to be the bilinear superpositions of the reduced amplitudes 
for producing a Ntt final state with allowed isospin by a charged isovector current: 



,AA' 



^AA' 
(7c = 



ttM |q|' 
2WQ^ 



( A\' , aX jA' \ 

(^0+^0+ + ^0-^0-) ■ 



The amplitudes for neutrino induced reactions are 

(7=3/2) 



(7=3/2) 
(7=3/2) 



(7=1/2) 
(7=1/2) 



Here x = ±3, ±1, 0± and only those resonances are allowed to interfere which have the same 
spin and orbital angular momentum as is in the following typical example describing the nTr"*" 
final state: 

^Ai (n7r+) Ai (n7r+) = 



E ai {Si) +V2j2ai (St,)] [E (^a^i) + ^^ E (^n)] 

+ E [E < (p^D + ^^ E «^ (Ay] [E (A"^) + ^ E (Ay 

i=i,3 

J=3,5 

+ E [E «i (^-35) + (FJ)] + ^/2 J] a^' (Pj); 



J=5,7 



^For the reader's convenience we use the same definitions and (almost) similar notation as in Ref. [5]. 
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Any amplitude (A^*) referring to one single resonance N* in a definite state of isospin, 
charge and helicity consists of two factors which describe the production and subsequent decay 
of the resonance: 

4 {N:) = {uN ^ N:) r){N: ^ Nn) ^ /^W r^W. 
The decay amplitudes, 77^'^, can be split into three factors. 



V 



sign(N:)Vx;^a(W), 



irrespective of isospin, charge or helicity of the resonance. Here, the first factor is the decay sign 
for resonance N* (see Table III of Ref. [5]), is the elasticity of the resonance taking care of 

the branching ratio into the nN final state and i]^l^{W) is the properly normalized Breit-Wigner 
term with the running width specified by the ttN partial wave from which the resonance arises 
(Eq. (2.31) in Ref. [5]). 

The resonance production amplitudes, f^^^\ can be calculated within the FKR quark model 
in exactly the same way as in Ref. [5]. It can be shown that they have the same structure as that 
given in Table II of Ref. [5] with the only important difference: the three coefficient functions 
S, B and C involved into the definitions of the amplitudes have to be modified. Indeed, since 
in our approach the structure of the polarization 4- vector e(\) has been changed with respect to 
that of the original RS model (by including the lepton mass and spin), we have to recalculate 
its inner products with the vector and axial hadronic currents. To do this, we used the explicit 
form for the FKR currents given by Ravndal [8]. As a result, the coefficients S, B and C (and 
thus the resonance production amplitudes) become parametricaly dependent of the lepton mass 
and helicity: 



B^B^ 



M 



6|q 



2 ' 



(Q(A)|q1--(A)-1U + 



aM 



Here 



(A) 



'-W- 



Q2 nQ Y 



aM ^ 3aM ) 



ZG^ {Q^ 
2W^|q*| 



W ' 2MW 

G^'^ (Q^) are the vector and axial transition form factors and the remaining notation is ex- 
plained in Ref. [5]. Other 5 coefficients listed in Eq. (3.11) of Ref. [5] are left unchanged. 
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Appendix: Coefficients y/% A\\ Kj\ 

Here we list the seventy nonzero coefficients \^/*^, A'^ and Kf^ appearing in Eqs. (8). 

3 ^ ^3^^ _ 2xfw^ + C[l + C^(l - 2xf]w'^ + 2(1 + C)xw, 
V^' = C\l - 2xfw' + [1 - C(2 - 0(1 - 2x)](l - 2x)w\ 
^35 ^ ^2(1 _ 23,)^3 ^ [1 _ - 0(1 - 23;)]«;2, 

14^^ = C(l - 2a;) W - (1 - C)(l - 2x)W, 
\/45 _ 2C(1 - 2x)w^ - (1 - C)(l - 4x)w^ 

= Cw' - (1 - 0^^^'; 

If ^ = 2(^X1^"^ + 2C(1 + C^)a;'«^, 

Vf 4 ^ 2(^X10^ + 2(1 - Cfxw, 

Vf 5 ^ 2(2(1 + 2x)xw'^ + 2(1 - Cfxw, 

^ 2Cx«;2 - 2(1 - C)xw, 
V^^ = ACxw^ - 4(1 - C)xw, 

V^^ = 4CVw^ + 2C[1 - 2C(1 - C)x]xw'^ - 2(1 - C)xw; 

Vf 3 ^ 20(1 _ x)w^ + 2C'^(1 - x)w - 1 - C, 
^ = 20(1 - x)w' + [1 - 2C(2 - 0(1 - x)]w, 
- 2CW - C(2 - 0^, 

vr = -2CV - (1 - o^ti^, 

= 2C(1 - - 2(1 - 0(1 - x)w, 
Vt = 2Cw' - 2(1 - Qw, 
Vt' = -2Cw' + 2{l-C)w, 

= C'w' - ('{1 - C)w', 

= 20(1 - 2x)w' - 2C [1 + C(l - 0(1 - 2x)] + 2(1 - 0^^, 
= 0(1 - 2a;)W - C [C(l - 0(1 - 2xy - 2x] - 2(1 - Qxw; 

Vi' = 2C'w' + 2Cil + e)w, 
Vt = 2CW + 2(1 - 0^^, 
K,35 _ 20(1 + 2x)w^ + 2(1 - Cfw, 
Vf6 _ _40x2ty2 - 2(1 - Cfxw, 
Vt = 2Cw' - 2(1 - 0^, 

vt = -^Cxw^ + A{l-C)xw, 

Vt = 4C^xw' + 2C[1 - 2C(1 - C):^;]^' - 2(1 - 0^^, 

Vt = 40(1 - 2x)xw^ - 4C[C(1 - 0(1 - 2x) + l]xw'^ + 4(1 - C)xw; 
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Af = C'(l - 2xfw^ + C[l + C'(l - 2xf]w^ - 2(1 - C)^;^, 

= C'(l - 2xfw' + [1 + C(2 + 0(1 - 2x)]{l - 2x)w\ 
Af = C'(l - 2x)w' + [1 + C(2 + 0(1 - 2x)]w, 
At" = C(l - 2xfw^ + (1 + C)(l - 2xfw^, 
Af - 2C(1 - 2x)w'^ + 2(1 + C)(l - 2x)w, 
Af = Cw + l + C; 

Af = 2C^xw^ + 2C(1 + C^)xw, 
Af = 2Cfxw'^ + 2(1 + C)^xw, 
Af = 2Cxw, 

Af = 2Cxw'^ + 2(1 + Qxw, 

Af = ew+e{i+o- 

Af = 2C^(1 - x)w^ + 2C^(1 -x)w + l-C, 

Af = 2C\1 - x)w' + [1 + 2C(2 + 0(1 - x)]w, 

Af^ 2(^10 + ({2 + 0, 

Af = -2CV - (1 + Cfw, 

Af ^2C{l-x)w' + 2{1 + (){!- x)w, 

Af = 2C^ + 2(1 + 0, 
Af = -2Cw' - 2(1 + O^', 

Af^C'w + Cil + C), 

Af = 2C=^(1 - 2x)w' + 2C[C(1 + 0(1 - 2x) - l]w - 2(1 + 0, 
Af = C'(l - 2x)^«;^ + C[C(1 + 0(1 - 2x)' + 2x]«;^ + 2(1 + 0^^;; 



4 33 
^5 


= 2CV + 2C(1 + C>, 


4 34 
^5 


= 2CW + 2(l + 0'«^, 


4 35 
^5 


= 2c^(l+xv+(l+o^ 


4 36 
^5 


= -4CV«;2 - 2(1 + O'a;^;, 


4 44 
^5 


= 2C«;' + 2(l + 0«', 


4 45 
^5 


= 2C«; + 2(1 + 0, 


4 46 
^5 


= -4:(xw'^ - 4(1 + C)xw, 


4 55 
^5 


= 2C^^ + 2C'(l + 0, 


4 56 
^5 


= 2C'(1 - 2x)w' + 2C{1 + 0(1 - 2x)«;; 
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-2C^(1 - 2x)W + 2C(2 - 3x)w, 
-C^(l - 2a;)W + 2(1 + C)(l - 2x)w, 
-C'(l-2x)w + 2(1 + 0, 
-C^(l - 2x) V + 2(1 - 0(1 - 2x)w, 
C(l - 2x) W, 
C{l-2x)w, 

-C\l - 2x)w'' - 2(1 - C)w, 
C(l - 2x)w', 
(w. 
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